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АВЗТВАСТ 

Design of linear control systems is examined using 
algebralc methods, Previous work indicates that the system 
describing equations are nonlinear, Various combinations of 
feedback and cascade compensated systems are analyzed 
using velocity constant, bandwidth and complex root speci- 
fications. Methods for linearizing;these system equations 
are demonstrated and examples are solved using specific 
digital computer programs, Feasibility of a general computer 


program is discussed, 
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I. з-с е е: [е 

Design is a word of many meanings and forms. Generally, 
it means a method or procedure for determining the value of 
system parameters such that the system will respond in a 
prescribed manner, These methods are usually derived for 
synthesis in the time domain by apolying modern state space 
techniques or in the complex frequency domain using root 
loes, Bade, Nyquist, or Mrerol"s plot. 

All the methods have their advantages and disadvantages. 
Time domain analysis with state soace allows the designer 
to examine the entire system at one time, Не сап only select 
one set of parameter values and investigate their effect on 
the system, The advent of digital computers greatly aided 
this process because the computer is easily adapted to a 
systematic trial and error design procedure, 

The complex frequency plane methods are inherently 
trial and error. These methods are graphical and extremely 
time consuming in their application, Present day programming 
is gradually coming to the rescue but it has limitations in 
the amount of useful and legible information that can be 
displayed on a graph. 

An outgrowth of the graphical methods was Mitrovic's 
method which was generalized by the Parameter Plane technique. 
This procedure transforms the analytical model of the system 
into a graphical representation. The benefit of this method 


is that it can handle multiparameters. 


Э 


From the parameter plane technique system design has 
returned to the analytical worid and developed a multi- 
parameter procedure called algebraic methods, Algebraic 
methods reduces the system differential equations and spec- 
ifications prescribing the desired response into a set of 
simu!taneous algebraic equations. The simplicity of the 
method is that it finds the parameter values that collect- 
ively provide the required system response, The method has 
developed over the years but its implementation on computers 
has not been explored, It is the intention of this thesis 


to investigate this area. 


II. Basic Theory 
The backbone of the algebraic methods is expressing the 
systems በ. чно equation as a polynomial of the form 
p 5 = © (2-1) 
ኗ= ር 


where s is the complex variable and is defined as 


5= > 5ሀ)ኳ + 3ሠሏህ]- 577 (2-2) 


and the ak's are real coefficients which are functions of 
the system parameters. Two independent equations in a,,S , 
and w, are obtained if equation (2-2) is substituted in 
equation ea) and the weal anid imaginary parts set edu 


to zero. Performing additional manipulation and introducing 


a new variable (1), (2), we obtain 


Y ew. f$ (9-0 (2-3) 


› ақық @x< (Š) =O (254) 

where 
9.5) з - Га5фл(5) + De-2(S)] (2-5) 

with 


6,9-0, 0Ц60--і ano Ø-A) (2-6) 
Other values of (Š) for various values of $ are tabulated 
in Appendix I, Observe that for a real root S= -@ and - 
substituting in equation (2-1) gives 


ነሰ K ሂ 
78 (1) 0,0-0 (2-7) 


— 


Previous work provided graphical methods for portraying 
root locations as a function of any two system parameters. 
(2), (3), (4), (5). If we are willing to abandon graphical 
analysis then the equations can include any number of unknown 
parameters. Proceed algebraically by substituting root 
locations in the basic equations (2-5) and (2-4) giving two 
equations as functions of all the unknowns, Requirement for 
solution is the number of equations equal the number of un- 
knowns. Cne obisii.g additional equations by one of two 
phi losophies; first, substitution of additional root locations 
into the basic equations; second, derivation of additional 
algebraic relationships from system performance specifications. 
Determination and solution of the proper number of equations 
results in a simultaneous set of parameter values that satisfy 
stated requirements. Naturally, if the parameter values are 
not satisfactory the performance specifications must be altered 
or the additional equations generated changed and a new 
solution affected, 

Inherently these equations provide a set of nonlinear 
simultaneous equations to be solved, The nonlinearities are 
generated from the system configuration or conversion of the 
performance specifications into algebraic relationships. Тһе 
purpose of this thesis is to develop methods of solving these 
nonlinear equations. Basic approach wili be that of linear- 
ization by substitution or manipulating the form of system 


components to facilitate linearization. When linearization 
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ІІІ, Third Order Development 

The basic system general block diagram is shown in 
Figure 3-1. In this section we will use a type one, third 
order plant to demonstrate the different solution techniques. 
Naturally, it is required that some type of compensation be 
employed, namely, feedback, cascade and feedback-cascade 
combination,so that the system satisfies the specifications. 
The third order system allows for development of the method 
while keeping the equation length reasonable, As needed 
reference will be made to higher order systems and some 
higher order examples will be solved, 
5.1 Feedback Compensation 

Consideration will first be given to tachometer feed- 
back, Figure 5-2, where the system gain (K) and the tachometer 
даіп (Ка) are unknown, Solution for the unknowns requires two 
simultaneous equations, Since the most common system speci- 
fication is a pair of dominant complex roots, equation (2-2) 
and (2-4) will be implemented, 


The system characteristic equation is 


c 
$°+(Р1+Рг]5 (PIP2* KK4)S * K=O (3-1) 
01:91109፡)111:1 (1:2 -ጩ.›፣፡ ከ... 1 ከ2. ео б апа Wr are 


substituted in the root equations giving 


кф ሟ፥ ቦጋዞ።ተ KK) A (S) (Річ Р) ф (5) + Шеф (8) =о (5-2) 


KOS) (PIP KKA Un A (PLA PAURAS (9-0 — G-» 


| 4 


It can be seen these equations are nonlinear due to the 
KK terme Linear equations are obtained by making the 


following transformation of variables 


X ak 
ү = КК Er 


After making the substitution and rearranging we get the 
following matrix form 

we ጊ 

А) Wa AO ]]x | [—piez2w. (4919204408) - мій фі (8) 

1 (5-5) 
` 2 ጋ 
ወር) ሠጪሊቆ፳|የ| |-ዞሀፌያሠለል(ዌ-ዐ4ሠ፻)ፀል(ኢ)-ኃ። 6.9) 
Linear equations are easily solved by diqital computers 
with programs using matrix algebra for implementinq the 
standard linear equation solving techniques. It is a simple 
Operation to inverse transform the linearizing substitution 
and obtain the parameter values 
K=xX 

and (5-6) 


Ks YI 


Observe that for this specific case the matrix equation 


degenerates to give 


yz 03685 Us) - x dz($ (827) 
Gina AG i 


y "n eL ee d G)-( PRPA Wn AlS- Pals) (3-8) 


P S) 


and 


because Ø, ($) = O 


Another possible method of feedback is acceleration 
feedback, Figure 2-3, which is very similar to tachometer 
feedback but affects the second order coefficient of the 
characteristic equation instead of the first order coeffic- 
ient, This type of compensation generates a system with two 
unknowns, therefore, requiring two equations which we will 
obtain from the root specifications. The system character- 


istic equation is 


S¥(PL+PZ+KKS + PLPZS+K 2 О (3-9) 


Appropriate substitution in the root equations yields 


KA S+PIP2 IM O. (S+ (Різ РенкХ, уш, @,\®9 2 ጴ(5) = ዕ (2208 


апа 


ኦሕጩ)ተየ[ሮረህነዩይን ፥!ክ!ቶ ዞደተ KKa) We ፊሼነተሠ። ቁጫ ። о (3-11) 


As in the tachometer feedback case the following linearizing 


transformation of vartables can be made 


Хек, 
VS ee 


Making the substitution and rearranging gives the matrix 


(3-12) 


simultaneous equations 


À (S) меф) ТХ) ሎጅ. (РаР Фе)” ихф(5) | (385 


ጺ(5) RARO Y| [PLA PEPA R QS) 
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For this type of feedback we do not qet the degeneracy 
duetothe^ortunate location of 2.05), but it is seen the 
matrix equations for tachometer and acceleration feedback 
are the same except for the second column of the coefficient 
matrix. Therefore, we can expect that when both types of 
feedback are used simultaneously the matrix equation will 
be very similiar to (3-5) and (3-12). At this point we will 
state the results if a fourth order B am is used, The 
matrix equation is 


PE) ሠ 065) | | K чу W. МУ урыннар ЧР Ba 
|| | - Мія" Фа| 5 


і 


መና መ ee 


m 


$43) DOS ке, E wn З) -(РІРо тРУРЗлРаРФун фо (РРО РУ фМЎ) 
ጋ 


ው wow dal) ( на 


From this it is seen that the form is the same as the third 


'፦ 


order case and all the increase in system order did was add 
a term to the riqht hand side, 

Let us now combine the types of feedback, Fiqure 3-4, 
and solve the problem. Immediately it is seen that before 
we can proceed we have to develop another equation besides 
the root equations. A frequently used specification is 
error constants which has different meaning depending on 
the type of system. For the type one system this means we 


can use the velocity constant defined as 


Қы зм በክ (5-15) 


with the steady state error for a ramo input, R(s) = I/s° 


Css = Lim SR(S) " Lim A. E | 


5-56. 1+ርፌ(5). 5-0 5645) ^ Kv (5-16) 
The characteristic equation for the system is 
CP. (PLePZ*KKa)S^- (PAP24KKX 4) S4 X20 Бы 
and the forward transfer function is 
боб) - a eae m (3-18) 
=[ ርፍትኮኒንርፍተኾይን ትኗጾቴ + КК] 
Apolying the root equations and the velocity coefficient 
equation we obtain the following set of simultaneous non- 
linear equations in K, KKy, and KK, 
KyPLP2 + Kv Ke -K =O (3-19) 
фа («К+ (PiP2+ ሾ ዩኒ ነህንሌ фь\$) “(атауы d (S) xU d. ነ = аф (5-20 


ፌን ኗ +ጺዮቷተ ኒን ሠፌፅ,ቢዔን +ሁኒተዞይትነ።ዩፌነጪዱቆፈኛን ሠሉ ቄኔ[81። ው (52) 


As in the separate type of feedback compensation we can apply 


the same transformation of variables to linearize the equations 
х= К 
Ро Ка 


Makino the substitution and rearranging we obtain the matrix 


equation | 
-4 бу б IK - КүРІР2 
фаз) YAE) MRS) JY] | -РІРашж dol) -CPL Pewa) -wihal | (9-22) 


i 


ኡ- 


| фіз) сім ФІЗ) АН Е F фас) (рата) фа) Waal) 


As predicted the above result is very similar to the 
feedback cases examined separately. The coefficient matrix 
second column is from the tachometer equation (3-5) and the 
acceleration equation (5-15). Similarly, the right side is 
not changed at all. From this we conclude the type of 
compensation for the root location specification affects terms 
in the coefficient matrix while the right hand side is affected 
only by the system order, 

A possible specification for generating an extra equation 
is bandwidth. The conventional definition of bandwidth will 
be used, that is, the frequency at which the magnitude of the 
output is one half the input magnitude for the closed loop 
transmission with S = j",. Because the magnitude involves 
squaring the transfer function we can expect to gain a 
nonlinear equation of a rather difficult form, For feedback 
compensation, it is found the problem is not extremely difficult. 

Consider the tachometer feedback, Figure 3-2, where the 


open loop transfer function is 





ЖӨН ШИТ 77 (3-2ዱ) 
(Рур) + KK. S 


and the closed loop transfer function is 


Gala oa JN (ЗОБ) 
SUSrPLGG+P2)+ KS + K 


to which we apply the bandwidth requirement giving 


Ke 


l- eee el T 
2." ፦ህህ LP1P2 -0 кке + [K-We (PLE PZ)} 
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To simplify the expression make the substitution 


A* -Ué (P2) 
and -21 
2 = Рірг – ше ЕР" 


whereby clearing gives 
—K +2AK - 2B KK tA WEB - To (3-28) 


Before proceding any further fet us examine the second equation 
required for solution. Select the velocity constant specifi- 
cation because a complex root specification generates two 
equations and it is not our intention to deal with over 


determined systems, hence using 


Ky = L. ҮҮ) = 
ፍጭር S(S+Pi)(S+P2) + KS 
gives 


Ky Kee K = —KyPLPZ (5-29) 


Examination of equations (3-28) and (3-29) shows an extreme 


nonlinearity. Fortunately, if we consider KK4 as an unknown 
we can use the elimination method to reduce the system to a 
quadratic equation in forward gain K, Solving equation (35-29) 


ት ነከ KK+ 


ረሄ 2 ዚ ደ КУР1Ў2)/ Ку (5-50) 


and substituting this value in equation (3-28) gives 


2 
[-\.0 -Wbh/k ie] t CIA -~ 2we/ey 4 2We PLPZ/KVIK 


(3-31) 
ж ді go pp АА ЕА y о зерви 2 
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This is easily solved for the system gain K yielding two 
values which can be substituted in equation (3-30) to give 


two values of Ky 


pire 
Kt, = hey p “д. (5-52) 


At this point the designer has a choice of which set of gains 
to use for the solution. Usually it is not a diff Т ОО се 
because one of the gain values is negative. Situations where 
both gains are positive values the designer can pick the values 
Beat со his pmmgeposiaem 

Мом we'!']] shift our attention to the acceleration feed- 
back, Figure 3-3. From previous expressions it can be expected 
the feedback gain in the second order term will help the 
solution, As in the tachometer case we will use bandwidth and 


velocity constant specifications. The velocity coefficient is 


ደሃ = E/prp2 
giving one equation 
Ka Ky PLPZ. (3-35) 
The bandwidth specification gives 
Ke 
fe — у —_ (3-34) 
2 ~ WAY Ді КК. т МВ 
ከ 

MEFE wAn pprc AND B= (РІ РА) 


Rearrangement and substitution for K from (3-33) gives @ 


a quadratic in Ка, 


21 


4 2 mg i % 2 
Wey (PLP2KV) Z. + (-ë ¿Je (PftP2 Kv) e ZU API P2 Kv ) Ko. 


^ [Bt eun 28/1 P2v – (РІ Р2 ки) °? = Ó (5-55) 
Aqain this qives two values for the acceleration and forward 
gains. As in the tachometer case the choice is left to the 
designer but usually one of the choices will have an unde- 
sirable negative gain. 

5.2 Cascade Compensation 

Cascade compensation is a widely used method of compen- 
sation, In general, it is commonly referred to as a lead, 
lag, or lead-lag filter, In this development we will not pre- 
specify the type of filter or a criteria for realizability. 
The decision of realizability is left to the system or filter 
designer. 

The first case examined is the single section compen- 
sator, Figure 5-5, In keeping with the major requirement we 
need three equations to solve for the unknown system gain (K), 
filter zero(Z), and the filter pole(P). If given a complex 
root specification and a velocity constant we can generate 
the required equations, A possible unknown is eliminated by 
considering any gain associated with the filter to be included 
in the overall system gain. 


The system characteristic equation is 


s*+ (Pls P24 P)S + (P(pa+P2) ppape)S+(APZPLK)S4KEZO — е) 


22 


indicating one nonlinearity, the KZ term. The velocity 


constant specification yields 


kv» . K(S*2) a hee — 
<+Р)(5+ғ2)(5+Р2) (5=0 Фір2Р 


which has the same nonlinear term as the characteristic 


equation, making the now familiar transformation of variables 
с (3-38) 


provides linearity with the three unknowns P, K, X which after 
affecting the root equations gives the matrix equation 


የጋ ~Ky PLP2 1 ы О 


| 
-FIP24 d, - (P4«-p2) 
қ ያ Фә й we da | 
2 2. 
wad, .- орған фу (ред owe fol |X|  — Perewdgy-(pt+P2) x 
|| . ሠ"ፏ-ሠፆፊ% (2-25) 


After solution we once again inverse transform the variable 


Wn Po PIP? Wa bot (PL+P2)Wa $) +l фа 24 Ж = 


г 
| 
| 
| 
| 
| 
| 
| 





for the filter pole 


Z ЮК (5-40) 


Unlike the feedback system a сһапде іп the system order 
affects not only the right side of the matrix equation but, 
also, the left hand coefficient matrix, The reason for this 
effect can be seen if the reduction of the system to a forward 
transfer function with unity feedback is examined, The filter 


pole multiplies all poles of the basic plant, thus inserting 
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menil ter pole in all the characteristic coefficients 
except a, and ag. 

The next logical consideration is a multisection cascaded 
filter. To demonstrate the best analysis form we will use the 
double section filter, Figure 3-6. The addition of each filter 
section adds two more unknowns, the section zero and pole, 
which requires two extra equations per section added, At the 
risk of completely specifying the system roots the easiest 
specification to use giving two equations is a root specifi- 
cation, otherwise other specifications will have to be trans- 
lated into algebraic equations. Тһе root specification can 
be very helpful since if not given explicitly it can be selected 
such that it does not effect dominance of any specified roots; 

Before proceding further, to ease the notation, multiply 


ቲከፀ plant denominator giving 


5 2 
where | 
| 
ል = Р, ት os 
É = Р, Pal (5-42) 
Reducing the overall system by normal methods, we achieve the 


following characteristic equation В 
= (А+? + Р2) 5% + ( Ву Я (Річ Р.) P1P2)5 * (B(P19 72) 


7 f 4 Ж ) 
+? (1 Pe) уе)» (ampaskdiokz2)S v Zella = O (22. 


At this point, note the transformation of the unknowns into 


the seven groupings PI, P2, PIP2, KZI, KZ2, К, апа К4122, TOE 


2^ 


double section may not be that unwieldy to solve but the 
situation gets worse as the sections increase, A simple 
manipulation can alleviate this undesirable condition, 
Conversion of the filter to polynomial form vice factored 
form orevents the generation of undesirable nonlinear terms, 
Granted some remain but they can be handled by appropriate 


variable transformation. The polynomial form of the filter 


15 
S^. XS 4 y 
S4 WS +Ë Eos. 
where X= 2,42 
T zi Z 
МЕ й № (5-45) 
Z= P 


Using this filter form in the system reduction we achieve a 


new and simpler characteristic equation 
- 2 
«З яз) (ВУ Qw 2)s”> (BWy A 2 +K ) s 


(2-26) 
(82 +0х)5 Ф КУ O 
The nonlinear terms can be linearized by 
R= KY 
-4 
жәй, (3-47) 


The application of the root equations twice and the velocity 


constant 


Kuss luam З. КС вк зау) ода (3-48) 
590 (520/52) (52,5485) Be 


£5 


yields the matrix simultaneous equation form 


< 
LF] E ss (3-49) 
ረ 
Wi 
where for root specification one the coefficients 

Аб) = dt) 
А(,234- А Wn À, (S) ша 909 4 Вал 46%) 
4(/, 53)  ሪህ4 #ሬ(ኛ) 
Alu 4)= aad (8) 

Aln S)= вам) з AU dal S) 4 Wide (9) 
AlZ i= Gol (3-50) 
д (2,2) = Вим Ф) 9и Д, (5) > ¿oO (S) 

4(5,3ህ e СӘ 

А 0,435 ша (5) 
А (2,5) = Вил“) АиФ) + 0“ фа (9) 

Bl) = ~ вои, (5) -Awt dals) є ам? tals) 

8%) = ~B Unh lS) -Awn Hl) -w dels) 


are functions of zeta one and omega опе, For root specifi- 


cation two the third 


rows one and two, ie, 


and fourth rows have the same form as 


alz i)z Alre) 


2 (4, 4) з 8(2,а) 


апа 


B(3) 
2/4) 


(3-51) 
= &)) 
8(2) 


М 
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where A(3,i), A(4,i), B(3) and B(4) are functions of zeta two 


and omega two, The remaining row is 


À (ZI) = -40 
А(52) т у В 
А(5с) = &(5) — ОО 2-3,4, 


(2-52) 


From this development we can see the polynomial! form 
filter provides the transformation of variables from Z; to 
Ка: , for Í = Í, ооо After solution of the linear equations 
the filter zeros and poles can be found by factoring the numer- 
ator and denominator. 

4.4 Combined Compensation 

Having considered feedback and cascade compensation 
individually we wifl now examine the combination of both types 
of compensation, Figure 3-1, Attention is given first to 
tachometer-cascade compensation, Figure 3-8, which indicates 
nonlinearity may be a problem. In this section the system 
gain is assumed known including the filter gain. Naturally, 
the three required equations are from velocity constant and 
root specification. 


The characteristic equation is 4 
З | rae 
stp (4+ P24 P)S* ሩ (PPL + ሥዖሥረ t PLP2 KI) S 


+ (APEP +KKE KISH KE =O (5-52) 


from which we apply the root equations and obtain the non- 


linear matrix equation 
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А5 с 2 E. O| 
| 
E cxi - L< (5-58) 
e О і ә 
E е | Ke | f 
where 
A= KLIP? 
B= «Ку 
— 5 


р= 2а фи) К 

£ = (22+ Р2) Wa DAS) tW bz 5) 

ሥ= ፖራ//ቹ (3-55) 
G = —PL?P2 ww dtt) -is PA) CC) - au hal) 
ሥ= tun’ dal!) K 

ድ = РА РАІ? «(PA рада фа б мо фа Ф) 

S = Ku, af 

T = -K n hi3) -PAP2U di 5)- Річ Рг) А 


Similar to the bandwidth situation this case lends itself to 
the elimination method of solution such that equation (3-54) 
reduces to a quadratic equation in tachometer gain K+ from 
which the compensator pole and zero may be found. As previ- 
ously the two answers т a choice which should operate under 
the same restrictions previously discussed, 

Finally, we shall consider the combination of accelera- 


tion feedback-cascade compensation, Experience indicates the 
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acceleration qain will 


simplify the equations, but for this 


case we find the compensator zero prevents this complete 


simplification. Using the root and velocity constant specifi- 


EU with the characteristic equation 
7, (P1 P2 4 Py Kn )3 9 (PAP2 VPP4A PP23 a B SŠ 
+ (PLP2ZP4K)S # KE =O 


(3-56) 
we obtain the nonlinear matrix equation 
re] d 
A 5 C ጋ Е ы 2 
D E | | (3-57) 
= E Р iii | 
A O = T [2 E UA 
where E т 
A= KPLP2 
B=-K 
C= (P14 P2) Wa pil) +e pals) 
= ¿a O 
= KU d.t) (3-58) 


{{ 


— P2 P2 (U^ dl f) = (PLHP) Wa pal S) -E dolf) 
- Kln hlf) = P1 Pe tua pel $) «(РРО «ША фа Р) 
a da($)-- 21 P2 D) v (Pt РІ) діл "фа (7) 
Аз э #7 
т = жол? (%) 


Тһе system of equations can be reduced by elimination (like 


li 


if 


H 


р 
o 
E 
G 
Ж 
ke 
= 


the tachometer case) to a quadratic in acceleration gain K 


from which the compensator pole and zero can be found 


In this section we see the simultaneous equations 


easily lend themselves to solution by elimination but appear 
to be close " Ebo ities 1S troe W^ a fourth S T uc 
is given that doesn't create nonlinearities (like bandwidth). 
The finearization method used throughout this thesis can be 
used for the K4Z or K4Z term, The situation is somewhat 
easier if a multisection filter is used, There is a nonlinear 
product for each section of filter, Therefore, if the number 
of Filter saati ons is even then specify an additional k/2 


complex roots, linearize and solve, 
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Third Order System With Tachometer Feedback 
Figure 3-2 
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Third Order System With Acceleration Feedback 
Figure 3-3 
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Third Order System With Tachometer And 
Acceleration Feedback 


Figure 3-4 
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Third Order System With Single Cascade Compensation 
Figure 5-5 
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Third Order System With Multisection Cascade Compensation 
Factored Form 


Figure 3-6 
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Third Order System With Cascade 
And Tachometer Feedback Compensation 


Figure 3-8 
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Third Order System With Cascade 
And Acceleration Feedback Compensation 


Figure 2-9 
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IV, Examples 

Examples in this section were solved on an IBM 360 
computer using subroutine programs furnished with the 
computer, The programs used are specific in nature and 
appear in Appendix II, The following points apply to all 
examples: 

a) negative feedback is desired so that negative 
feedback constants are not acceptable since this implies 
positive feedback. 

b) when a cascade pole or zero is found positive 
sign indicates left half plane and negative sign indicates 
right half plane, 

с) for bandwidth and feedback-cascade combination 
examples, there are no logic statements in programs to select 
the proper set of parameter values for determining the system 
roots. This part of the program was written with a priori 
knowledge. For roots the sign convention mentioned above does 


not apply, use the standard sign convention, 
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Example 4,1 


R V За ረመ Кк ር 
P ж. (S45) (S46) 
«ተፀ 
Figure 2,1 


Using feedback compensation design a system with 
open loop transfer function in Figure 2,1. 
Root Specification: Zeta = 0.5, Natural Frequency = 10.0 


No Velocity Constant or Bandwidth Specification 
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INPUT DATA 


kk aR ak i ak ok 

POLES ARE ZERO, 0,.50000D O1., 0.50000П 01 
КОСТ SPECIFICATIONS ARE ZETAs 0.500000 00, 
NATURAL FREQUENCY= 0.10000П 02 


NO VELOCITY CCNSTANT SPECIFIFD 


NO BANDWIDTH SPECIFIED 


QUTPUT DATA 
Ж жж жжҗжк Жж 


FORWARD GAIN 


О. 10000Е 02 


TACHOMETER GAIN 


С.80000Е OC 


HARACTERISTIC EQUATION 


-о 


0.10000Е 03 9.11С00Е 03 0.11000Ғ 02 


THE SYSTEM ROCTS ARE 


RFAL PART 


-0.10000E 01 -0.50000E 01 -0.50000E 01 


IMAGINARY PART 


0.0 -0.86603Е 01 0.86603F 01 


De 


0.10000Ғ 01 


Examole 4,2 





Figure 4,2 


Design an acceleration feedback system with the open 
loop transfer function shown in Figure 4,2, 
Root specification: Zeta = 0.5, natura! frequency = 1120 
No velocity constant and bandwidth specification. 
This examole indicates an unattainable solution since 
solution gives negative K and Ka anda ЕИ half plane 


moO. 
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ІМРОТ ОАТА 


OIITITIITII: 

POLES ARE ZERO; 0.100000 01, 0.500000 01, 
ROOT SPECIFICATIONS ARE ZETA= 

NATURAL FREQUENCY= 0.110000 


NO VELOCITY CCNSTANT SPECIFIED 
NO BANDWIDTH SPECIFIED 


OUTPUT DATA 
ROR I RR ak ak 


FORWARD GAIN 


->0,100105 ዐፊ 


ACCELERATION GAIN 


-0.82645E-01 


HARACTERISTIC EQUATION 


ч- СУ 


THE SYSTEM ROCTS ARE 


REAL PART 


-0.34194E 01 0.24194Е 01 -0.55000F 01 


IMAGINARY PART 


0.0 0.0 -0.95263E 01 
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0.500000 00, 
02 


1 0.30000Е 02 0.123T3E 03 


0.600000 01 


0.12000Е 02 


-0.55000Е 01 


0.95263E Ol 


Example 74.5 
-Қа- (2-тен Ff 
, $(S+1)(S+5) (S+6) 


Kas? 


Figure 4.3 


Redesign of Example 4.2. 
Root Specifications: Zeta = 0.5, Natural Frequency = [2.0 


No Velocity Constant and Bandwidth Specified 


4 | 


INPUT DATA 


ተተተተተቸቸፋችቸችት 

POLES ARE ZERO, 0.100000 01, 0.500000 01, 
ROOT SPECIFICATIONS ARE ZETA= 0.500000 00, 
NATURAL FREQUENCY= 0.120000 02 


NO VELOCITY CONSTANT SPECIFIED 
NO BANDWIDTH SPECIFIED 


OUTPUT DATA 
жжжжжжжжжжж 


FORWARD GAIN 


0.36000Е 03 


ACCELERATION GAIN 


0.29306Е СС 

ር ው EQUATION 

i 0.30000E 02 0.14650Е 03 
THE SYSTEM ROCTS ARE 


REAL PART 


-0.60000E 01 -0.60000E 01 -0.37068E-07 


IMAGINARY PART 


-0.10392Е 02 0.10392Е 02 -0.15811Е 01 
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0.600000 01 


0.12000Е 02 


-0.37068Е-07 


0158 ПЕ 01 


Examole 4.4 





Figure 4,4 


Design a third order system with acceleration and 
tachometer feedback for the open loop system in Figure А 
Коої Specification: Zeta = 0.7, Natural Frequency = 10.0 
Velocity Constant = 6.0 


No Bandwidth Specification 


45 


POLES ARE ZERO, 0.100000 01, 0.200000 01 


0.700000 00; 


VELOCITY CONSTANT IS 0.600000 01 
NO BANDWIDTH SPECIFIED 


OUTPUT DATA 
ak a ak ak ake k KR 3k OR 3R OK 


FORWARD GAIN 


0.37500Е С4 


TACHOMETER GAIN 


0.16613Е OC 


ACCELERATION GAIN 


0.12933Е-01 
COMPUTED VELOCITY СОМ5ТАМТ IS = 0.60000Е 01 
олады EQUATION, 
0.37500Е 04 0.62500Е 03 0.51500Е 02 


THE SYSTEM ROCTS ARE 
REAL PART 


-0.70000Е 01 -0.70000Е 01 -0.37500Е 02 
IMAGINARY PART 


05 715155 0]1 -0.7101%Е Ol 0.0 
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0.10000Е 01 


Ехатріе 2,5 






4  — Ko ፡፡ - 
9(8-2)(8-6)(5-0)6-10 







KaS“+K+S 





Figure 4.5 


Design a fifth order system with acceleration and 
tachometer feedback for the open !oop system Figure 1.5 
Root Specification: Zeta = 0.5, Natural Frequency = 5.0 
Velliecity Constant = O85 


No Bandwidth Specification 
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INPUT DATA 
222213424; 


POLES ARE ZERO, d 
0.500000 00, 
01 


VELOCITY CONSTANT IS 0.500000 00 
NO BANDWIDTH SPECIFIED 


OUTPUT DATA 
X kc 3k OK 


FORWARD GAIN 


0.13349Е 04 


TACHOMETER GAIN 


0.11909Е 01 


ACCELERATION GAIN 


0.65497Е-01 


COMPUTED VELOCITY CONSTANT IS = 0.50000E 00 
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ТНЕ 5Ү5ТЕМ КОСТ5 АКЕ 


REAL PART 


ዛ -0.11184Е 02 
IMAGINARY PART 


О. 34641Е 01 
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ч- (^) 
me 
P 
2) 
> 
С 
ጫዳ 
m 
х 
ፍመ 
ዕ/ነ 
መ 
- 
СУ) 
m 
E 
> 
- 
= 
= 
2. 


0.19354Е 04 0.25400F 03 


-0.11184Е 02 -0920000E 01 


0.26026Е 01 -0.3464176 07 


Example 4,6 





Figure 4,6 


Design a tachometer feecback system for the open 
looo system in Figure 4.6., 
Мо Root Specification 
Velocmuy seencrant = 0.9! 
Bandwidth Frequency = 4.0 
Systems roots are for number two gain values since system 


gain one is negative, 
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POLES ARE ZERO, 0.50000Е 01, 0. 60000Е 01 


NO ROOT SPECIFICATION 


VELOCITY CONSTANT [፲ 5= 0.91000F 00 
BANDWIDTH FREQUENCY І5- 0.40000F 01 
OUTPUT DATA 
TIITTTPTIT 
SYSTEM GAIN ONE= -0.31554E 02 
TACHOMETER GAIN ዐርክር= .20496Е 01 
SYSTEM GAIN TwC= 0.,41920E 02 
TACHOMETER GAIN TWO= 0.38325Е 00 
COEFFICIENTS OF CHARACTERISTIC EQUATION 
(ASCENDING ORDER) 

0-41920E 02 0.46066Е 02 0.11000Е 02 0.10000Е 01 


THE SYSTEM ROOTS ARE 


REAL PART 


-0+12317TEn01 -0.48841Е 01 -0.48841E Ol 


[IMAGINARY PART 


0.0 -0.31905Е 01 0.31905Е 01 
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Example 24.7 





` —— 
—— 


ГҮСТЕ | 5127)“ 


Figure ЖТ 


Design an accelerometer feedback system for the Open 
loop system Figure 4.7. 
No Root Specification 
Velocity Constant є 5,0 
Bandwidth Frequency = 2,0 | 


System roots are for gain values number one since Ka2 


is negative. 
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POLES ARE ZERO, 0.30000Е 01, 0.40000Е 01 
NO ROOT SPECIFICATION 
VELOCITY CONSTANT IS= 0.50000Е 01 


BANDWIDTH FREQUENCY І5- 0.30000Е 01 


= 0.60000Е 02 
IN TwOz -0.50494Е 00 


HARACTERISTIC EQUATION 


መርን 


0.60000Е 02 0.12000F 02 0.99633E 01 0.10000Е 01 


THE SYSTEM ROOTS ARE 


REAL PART 


-60.29863E OC -0.2986ЗЕ 00 =0„9З6б60Е OI 


IMAGINARY PART 


0.25134Е 01 -0.25134Е 01 0.0 
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Example 4,8 











К 
(SPs 





Figure 4,8 


Compensate the open loop system in Figure 4.8 using 
a single section cascade compensator, 


Root Specification: Zeta = 0.5, Natural Frequency = 6.0 
Velocity Constant = 4.15 


No Bandwidth Specification 


5 


POLES ARE ZERO, 0.500000 00, 0. 700000 01 
ር ZETA= 0.500000 00, 
0. 600000 01 


VELOCITY CONSTANT IS 0.415000 01 


NO BANDWIDTH SPECIFIED 


OUTPUT DATA 
жжжжжжжжжжж 


SYSTEM GAIN 
0.13756E 04 

COMPENSATOR PDLE 
0.37847Е 02 


COMPENSATOR ZERO 


0.39961E OC 
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VELOCITY CONSTANT(COMPUTED) 


0.41500E 01 


COEFFICIENTS OF CHARACTERISTIC EQUATION 
(ASCENDING ORDER) 
0.54972Е 03 0.15081Ғ 04 0.28735Е 03 
0.45347Е 02 0.10000Е 01 


THE SYSTEM ROCTS ARE 


REAL PART 


-0.39200Е CC -0.30000Е 01 -0.30000E 01 


IMAGINARY PART 


0.0 -0.51962Е 01 0.51962Е 01 
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-0.38955Е ШЕ 


INPUT DATA 

* doe d oe s x é$ k 

POLES ARE ZERO; 0.100000 01, 0.150000 02 
ROOT SPECIFICATIONS ARE ZETA= 0.665000 00, 
NATURAL FREQUENCY= 0.596000 01 


VELOCITY CONSTANT IS= 0.280000 02 
NO BANOWIDTH SPECIFIED 


0.506000 00, 
02 


FORWARD GAIN 


0.3463ТЕ 04 
FILTER ZEROS ARE 


REAL PART 
-0.51508E 01 -0.91788Е 01 


IMAGINARY PART 
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Example 4,9 









Á Чы 
sms tS) 





а нн away ፥ጻ ባኩ 
aaa puna Айы уна чр чего LO y ጭሌ ж, ыр Aaa ect. ING? es ee ee eee 


Figure T» 


Design a double section compensator for the opok 
loop system Figure 4,9, 
Root Specifications: Zeta | - 0,665, Natural Frequem MEM 
5,96; Zeta 2 = 0.506, Natural Frequency 2 = 14.66 
Velocity Constant = 28.0 


No Bandwidth Srecification 


FILTER POLES ARE 


REAL PART 


-0.14107E 02 


IMAGINARY PART 


-0.1381ТЕ 02 0.1381ТЕ 02 


HARACTERISTIC 


መርን 


THE SYSTEM ROCTS ARE 


REAL PART 


-0.39634Е 01 


| | 
eo 
ее 

O 
ዮ“ ም 
оо» 
mm 
Oo 
ht 


[IMAGINARY PART 


0.44512Е 01 
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EQUATION 


0.1 
1 


910123 ዮ 05 
9094100005 01 


-0.21451Ғ 02 -0.74180Е 01 


0.0 -0.12645Е 02 


Example 4.10 


21010 
з (5-5 ) (5-8) 





Figure 2,10 


Design a single section cascade compensator with 
tachometer feedback for open loop system Figure 5.10, 
Root Specification: Zeta = 0.5, Natural Frequency = 7.0 
Velocity Constant = 0.5 
No Bandwidth Specified 
These specifications can not be satisfied, Choice one 
gives negative feedback gain and roots indicated, Choice 


two gives undesirable right half plane compensator. 
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0.80000Е 01 


0.50000Е 01, 


SYSTEM POLES ARE ZERO о 


0.50000E 00, 


E ZETA= 
0. 10000Е 01 


0.50000Е 00 


VELOCITY CONSTANT IS 


NO BANDWIDTH SPECIFIED 


CHOICE ONE 


CHOICE TWO 


TS OF CHARACTERISTIC EQUATION 


0.20664Е 03 0.28640E 02 


0.11035Е 04 


THE SYSTEM ROOTS ARE 


REAL PART 


-0.35000E 01 -0.35000Е 01 30 21352Е 02 


-0.28840Е 0С 


IMAGINARY PART 


0.0 


0.60622Е 01 


0.0 
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Example 4,11 


ego 










200 
5(8+5) (8+8 


KTS 


Figure 4,11 


Redesign Example 4,10 for single section cascade 
compensator with tachometer feedback for open loop system 
Figure 4.11. 

Root Specifications: Zeta = 0.5, Natural Frequency = 10.0 
Velocity Constant = 0.5 

No Bandwidth Specified 

Data on page 6! shows results for Wn = 10,0 and system 
roots for parameter vaiues choice one. [f Wn is changed 
to Wn = 14,0 data on page 62 shows both choices of 


parameter values are satisfactory. 


60 


NPUT DATA 
ተቹት ፍች 3k a K 


I 
* 


0.80000Е 01 


0.50000Е Ol, 


SYSTEM POLES ARE ያርጸዐዩ፥ 


0.50000Е 09, 


ZETA= 
10000Е 02 


Е 
0. 


0. 50000Е 00 


VELOCITY CONSTANT IS 


NO BANDWIDTH SPECIFIED 


ONE 


CHOICE 


CHOICE TWO 


F CHARACTERISTIC EQUATION 


0.16483F 03 0.16409E O2 


0.64829Е 03 


THE SYSTEM ROCTS ARE 


REAL PART 


-0.62909E Ol -0.50000Е 01 -0.50000Е 01 


-0.11786Е 00 


IMAGINARY PART 


0.86603E 01 


-0.86603E 01 


0.0 


0.0 


6 | 


0.80000Е 01 


0.50000Е 01, 


SYSTEM POLES ARE ZERO,, 


0.50000Е 00, 


Ах 
00Е 02 


= 0.50000E 00 


VELOCITY CONSTANT IS 
NO BANDWIDTH SPECIFIED 


CHOICE ONE 


CHOICE TWO 


TS OF CHARACTERISTIC EQUATION 


N 
G ORDER) 


We 


0.24Т18Е 03 0.17616Е 02 


0.71652Е 03 


THE SYSTEM ROOTS ARE 


REAL PART 


-0.34564Е 01 -0.70001Е 01 -0.70001Е 01 


-0.15980Е 00 


IMAGINARY PART 


0.12124Е 02 


-0.12124Е 02 


0.0 


0.0 
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Example 4,12 


p UU И 
5 (8+5) (8+8) 





Design single section cascade compensator with 
acceleration feedback for open loop sy-tem Figure 4.12, 
Root Specifications: Zeta = 0.5, Natural Frequency = 5.0 
Velocity Constant = О а 
No Bandwidth Specified 


System roots are for choice two parameters 
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0.80000Е 01 


0.50000Е 01, 


SYSTEM POLES ARE ZERO,, 


0.50000Е 00, 


As 
00Е 01 


0.50000Е 00 


VELOCITY CONSTANT IS = 


NO BANDWIDTH SPECIFIED 


CHOICE ONE 


CHOICE TWO 


EQUATION 


CHARACTERISTIC 


0.24345E 03 0.46322Е 02 


0.10923Е 0% 


THE SYSTEM ROOTS ARE 


REAL PART 


-0.25000Е 01 -0.25000Е 01 -0.41033Е 02 


-0.28868Е 00 


IMAGINARY PART 


0.0 


0.43301Е 01 


-0.43301Е 01 


0.0 
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V. General Proaram 

The orobiems salved Rm this thesis were donc orhan 
individual basis with no attempt made to write a genera! 
program. In retrosoect it is seen that a general orogram 
couldabe written oprowided the typeof specifica On R 
ре satisfied are selected, the maximum system order Inciudina 
any -comoensation effects is stated, the desired output 
information and format is known, and the system configuration 
is consi dered, 

Qur-prime specificatWon"was-^a comstex root spec!|fwWwe 
cation which is dependent оп thes form of the characteris Tie 
equation coefficients, In all the cases consi dered the right 
hand side of the matrix equation had the same form, The form 
is a function of system tyne and order, and has two general 


morms 
V, Г 


com. ኒ o 
Z. - (k QA be. (S) 


mem 


кп (5-!) 


from equation (2-5) апа 
A 


7 as iR 4206) 
б (5-2) 


from equation (2-4) where 
em О 


ድ = a PA oni Pio Vs: e 


z ኣ 
5 


(ри-с сым ә bass sor б Ман . ‚ (5-3%) 


Nee Vert hand side coeficient matrix aiso reduces wore 


convenient form 
= J КЕ =. 
| 7199 ша 4%(%)---- LU $,., 09) №, | 
i |. ንና | 
ያረ( е) wr d if) FE (UA фу 69 |ኖ/ | 


ቲክ coefficient 


where R is an unknown appearing in the j 
of the characteristic equation. If a particular coefficient 
of the characteristic equation has no unknowns then the 
appropriate place in the coefficient matrix above is replaced 
by a zero. As noticed for the cascade comoensation system 
terms similar to the right hand side are aenerated in the 


left hand coefficient matrix multiolyina the unknown filter 


pole or polynomial coefficient except the coefficients are 


п ኑ 
е On Na ф, (9) 


кох (5-5) 
апа 
ነ 
< 
ሷኗሠዚ Ф.. (8) 
K- 
5 | (5-6) 
sre K= n+ ORDER OF FILTER 
ге ^o 
= 
n-i = ae p 
Ди. = Sum or THE PRopucr OF ¢ OF THE FILTER Poles 


¿ = /, Ре K-n 
The root equations can be generalized by applying them 


to a characteristic equation with coefficients 


Ак = к х + bk É + CK (5-7) 


and then the result can be reduced to a given situation by 
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reading in the orooer value for the coefficients &« , 6, and cr. 
The velocity constant also has certain properties that 
таке it conducive to programming, The basic form is 
ዕ 
Ky E ረ ШЕ с 
፪ - І (5-3) 
// ќа 
where n is the number of system poles and m is the number of 
system zeros. Тһе form above aoplies for systems with unity 


feedback and acceleration feedback. If tachometer feedback 


is used the velocity constant becomes 


m | 
Ios m KO 


(5-0) 


з ч жен e aes o EA 





i w 
J P+ KAQ 77 Ze 
The remaining situations are for cascade? systems an^ 
the form is now 
በሰሞ 
Ky = K WL / ሯ/ 
ra. x 
ГА 4. 


where іФ їПпе пипбег ої ‘cascade filter sections. The cascade 


(5-10) 


acceleration combination is as above but the tachometer 


cascade combination becomes 


menm A. 
Ky = Kite Mes AA... (5-11) 





In al! the above situations if any pole or zero is not in the 
Зетет іі is reolaced by a one where as removal of tachometer 
feedback imolies K. = nor 

The final specification presents a rather tedious oroblem 


as the general form shows 
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J| Gina S 7% ----- у йо) 


—— = 


ነ- “ሸህ. ИСА е И 
> (an-ı Sn An-3S+---+1)+ (dnS+an-28 + .-4%) $24, 5-7 12) 





І 
= E 
where the a,'s are a function of the unknowns. Since the terms 
are squared it is recommended that bandwidth be set up for the 
most difficult case exnected to be encountered. 

Once the soecifications have been oreset in the program 
they can be implemented by reading in aporopriate constants, 
then using logic statements the simultaneous matrix equations 
can be formec, Again as in the bandwidth case the nonlinear 
cases using the elimination method for solution should only be 
set up for the maximum system configuation and order expected 
to be encountered. This is because the computer can process 
numbers but is unable to process unknowns specified by alpha- 
betic characters. Designers may find that for cases of this 
nature it may be beneficial to use an aporopriate iteration 
routine mentioned in the next section, A rough general flow 


chart of major operations is shown in Figure 5-1. 
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Read Data 





ጨፌ» -- E ЕСШЕ ЕЕРЕЕ 


И 
IPrir* Input Datal 
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Specifications 







Velocity 
Constant 






Bandwidth 


| 
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Equations 






EE ES i ቁ “ቁ መም: ማመ ^o AD, Idi AL Rn OM nawa wa —rmsarsaaarmÑrCtnpapasamaaapaas — a, 
= bog ye з с Ee LR TT Re^ ap чаас. гатары ate сор 


Еос Та опао г (ጩህ. 
Speci ፡ቭኛ ት Cheek 
Other Desired Computatio 


eet aa... ЧА ЬО О desc co am 


ብ 





лылық. ወዬ «т. 


Print Output 


General Program F lowagraph 
Figure 5= | 
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VI. Conclusions and Recommendations 

The purpose of this paper was to investigate the feasi- 
bility of using a computer for designing with algebraic methods, 
It has been shown that this can be accomplished, in most cases, 
by the linearizing transformation of variables or manipula- 
tion of system analytic form, This was not true for the 
feedback-cascade compensation cases or when bancwidth speci- 
fication was to be satisfied but these oroblems could be 
solved by application of a straight elimination process. 

The important advantage of this method is it only restrains 
the system sufficiently to solve for the unknowns, Other 
design attempts tended to specify all root locations in imple- 


menting a solution, 
An area of possible future study is converting other 
system specifications into algebraic equations and analyzing 


techniques of solution required. We have dealt with only 
three specifications in this paper. Can we effectively use 
peak overshoot, steady state error, settling time, resonance 
frequency, or rise time? In analyzing these specifications 
nonlinearities may aovear that defy solution by elimination, 
therefore, it may require, u5simgs іже певні on methods for 
solution like Newton*s method or steepest decent, or a 
combination of these. The results obtained by these methods 
are dependent heavily on the initial value guess and to a 


lesser degree the quality of convergence to the final value, 
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An interesting ооіпті ТОҒ Тигіпег іпмесітсегітсот әш С 
be the possibility of deriving another definition for band 
width based on a criteria other than magnitude. If band- 
width can be specified as a comolex number then the non- 
finearities are not compounded by squaring and the specifií- 
cation wil! reduce to a two equations similar to the root 


equations. 
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